A Re-examination of Quenches in 4 He (and 3 He) 
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In the light of recent difficulties in observing vortices in quenches of liquid 4 He to its superfluid 
state we re-examine the Zurek scenario for their production. We argue that experiments in He are 
unlikely to produce true vortices in the numbers originally anticipated, if at all, because of the wide 
Ginzberg regime and the slowness of the mechanical quenches. On the other hand, the observed 
production of unambiguous vortices in neutron-bombarded 3 He, with its narrow Ginzberg regime, 
and fast quenches, is to be expected. 

PACS Numbers : 11.27. +d. 05.70.Fh, ll.10.Wx, 67.40.Vs 



In thermal equilibrium the behaviour of simple systems 
experiencing a continuous phase transition is generic, as 
manifest in the utility of Landau-Ginzburg theory. Is this 
also true dynamically? The relevance of this question is 
that the early universe proceeded through a sequence of 
phase transitions whose consequences are indirectly ob- 
servable, but whose detailed dynamics is unknown. Al- 
though it is difficult to measure an order parameter as 
it changes, many transitions generate topological charge 
or topological defects which can be detected, and from 
which the evolution of the field can be inferred. Moti- 
vated in part by Kibble's mechanism for the formation of 
cosmic strings in the early universe , Zurek suggested 
H that we measure the density of vortices produced dur- 
ing quenches of liquid A He and 3 He into their superfluid 
states. 

The original scenario, as proposed by Zurek, is very 
simple. It is exemplified by assuming (as is roughly 
appropriate for He) that the dynamics of the transi- 
tion can be derived from an explicitly time-dependent 
Landau-Ginzburg free energy of the form 



F(t) = / d 6 x 



2m 



*(W| 2 + i/3H 
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In ([j]) <p — (4>i + i4>2)/ \/2 is the complex order-parameter 
field, whose magnitude determines the superfluid den- 
sity. In equilibrium at temperature T, in a mean field ap- 
proximation, the chemical potential a(T) takes the form 
a(T) = a>oe(T c ), where e = (T/T c — 1) measures the crit- 
ical temperature T c relative to T. In a quench in which 
T c or T changes we assume that e can be written as 

e(t) = e - — 9{t) 

for — oo < t < tq(1 + eo), after which e(t) = — 1. eo 
measures the original relative temperature and tq defines 
the quench rate. The quench begins at time t — and 
the transition from the normal to the superfluid phase 
begins at time t = cqtq. 

With £q = h 2 /2mao and To = h/ctQ setting the funda- 
mental distance and time scales, the equilibrium correla- 



tion length £ e ,j(Ai) and the relaxation time r(At) diverge 
when the relative time At = t — cqtq vanishes as 



£ eq (At) - Co 



At 



-1/2 



r(At) = r 



At 



Meanwhile, the non-equilibrium correlation length £(At) 
at relative time At is defined from the diagonal equal- 
time correlation function (up to powers) 

(Mr)M0))M = S ab G(r,At) ~ <Ue" r/e(At) , 

for r large compared to £ (but not necessarily asymptot- 
ically large). As we approach the transition, eventually 
the relaxation time will be so long that the system will 
not be able to keep up with the pressure or tempera- 
ture change. We estimate the relative time —t < at 
which the change from equilibrium to non-equilibrium 
behaviour occurs by identifying t(— t) with t i.e. t = 
y/TQTQ. After this time it is assumed that the relaxation 
time is so long that £(At) is more or less frozen in at 
the value £_ := f(-t) « £ e9 (-t) = £o(tq/t ) 1/4 until the 
system is again changing slowly, at time At ps +t, when 
it attains the value £+ := £(+t) ~ S, eq {+t) = £,eq{—t)- 

If, for the sake of argument, we assume that the field 
fluctuations are approximately Gaussian, then the field 
phases are correlated as (e ie ^ e- ie ^) A t ~ e - r /^ At \ 
where o (r) = |0 a (r)|e i9 W, on the same scale as the 
fields. The correlation length of the field can only be 
measured indirectly, but the field fluctuations will en- 
code vortices as defects in the phase. If the initial vortex 
separation (in a plane) is £de/, then their density, ridef, is 
n def — 0(l/^ e f). Zurek initially made the assumption 
that £def ~ £+ whereby, with £ + « in turn, 



ndef 




(2) 



Since £o also measures cold vortex thickness, tq 3> to 
corresponds to a measurably large number of widely sep- 
arated vortices. We note that, since the argument is ba- 
sically one of causality, a free energy of the form ([[]) is 
not strictly required. 
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The simple order of magnitude prediction (0), with RG 
improvement for He, has been the motive for several ex- 
periments. In vortex production in quenches of 3 He into 
its £?-phase two separate experiments at Grenoble j| and 
Helsinki show quantitative agreement with (|J) and one 
experiment || at Lancaster on vortex production in 4 He 
shows consistency. However, a second experiment (tJ on 
A He finds no vortices. Because of this, it is timely to re- 
examine this picture, with its implicit assumption of an 
initial domain structure characterised by a single length. 
In this brief note we shall argue that the prediction (|^) 
is unreliable because vortex separation cannot be related 
simply to field or phase correlation length. However, the 
success of the 3 He experiments shows that there can be 
a de facto similarity, and we shall show how this arises. 

To see this we need an explicit model, and again we 
adopt (0). Motivated by Zurek's later numerical ||] 
simulations with the time-dependent Landau-Ginzburg 
(TDLG) equation for F of ([j]), we assume a linear re- 
sponse 



r dt 



SF 



Va, 



(3) 



This does not lead to an asymptotic fall-off of the form 
e- r ' D . In fact, for large r, G(r) cx exp(-0((r/L>) 4 / 3 )). 
Nonetheless, numerically, it is remarkably well repre- 
sented by e~ r l D , with coefficient unity in the exponent, 
for r being a few multiples of D, for reasons that are not 
clear to us. In that sense Zurek's prediction for a cor- 
relation length of the form D — ^(tq/tq) 1 ^ 4 is robust, 
since explicit calculation || shows that does not vary 
strongly in the interval — i < At < i. 

Despite this approximate constancy, the relationship 
between the unobservable and the observable ^ef 
is, a priori, very weak since the two measure very differ- 
ent attributes of the field fluctuations. Whereas the core 
of every vortex is a line zero of the complex field <j> the 
converse is not true since zeroes occur on all scales. How- 
ever, a starting-point for counting vortices in superfluids 
is to count line zeroes of an appropriately coarse-grained 
field, in which structure on a scale smaller than the 
classical vortex size, is not present fTo(| . This is, indeed, 
the basis of the numerous numerical simulations |ll]] of 
vortex networks built from Gaussian fluctuations (but see 
p"2[). For the moment, we put in a cutoff I = O(£o) by 
hand, as 



where rj a is Gaussian noise. We go further and assume 
that, for early times, prior to the completion of the sym- 
metry breaking at least, the self- interaction term can be 
neglected (f3 = 0). This both preserves Gaussian field 
fluctuations and leads to £± « ^{tq/tq) 1 / 4 arising in a 
natural way, as we shall see. 

Firstly, we anticipate that correlation lengths are ap- 
proximately frozen in during the time interval —i < At < 
i. Thus, in the first instance it is sufficient to perform 
our calculations at At = 0, for which D := £(At = 0) « 
w Writing the resultant equation in time and 
space units tq and £o as 



G(r) 



fke lkx G(k)e 



We stress that the long-distance correlation length D de- 
pends only on the position of the nearest singularity of 
G(k) in the complex k-plane, independent of I. 

This is not the case for the line-zero density n zero , 
depending, in our Gaussian approximation jl3 14 1 , on the 
short- distance behaviour of G(r), 
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-1 G"(0) 

2ne zero ~ 2vr G(0) : 



(5) 



0«(M) = - 



6 - —9{t) 



a (k, t) +T T] a (k,t) 



the ratio of fourth to second moments of G(k) e k 1 . In 
evaluating 



gives the solution, at At = 0, a (k) 



tq I dt exp 



dt' 



e o - — 0(t) 



%(k). 



The resulting un-normalised correlation function has 
power (Fourier transform) 



G(k) = I dtexp 



dt' 



k 



e - —0(t') 



(4) 



For a typical quench in 4 He, e$ ~ 10~ 2 is very small, but 
tq ~ 10 10 is so large that e$TQ, e^tq » 1. G(k) is then 
approximately independent of e , 

fT Q (e +k 2 ) 

G(fc) = e T Q fc / dte~ t/T °. 



dk k 2 e 



G(k) 



dkk 4 e' k l ' 



G(k) 



we substitute for G(fc) from (Q) and perform the k inte- 
gration first. For small eo, very large tq, and I = 0(1), 
in dimensionless units the dominant contribution is from 
t ~ £otq. On neglecting terms relatively 0(e~ e ° TQ ) we 
find, at At = 0, 



e -t /tq I r°° 

x / di FTWP/l dt 



[t + P/2] 5 / 2 



0(l 2 ) 
(6) 



independent of eo- The details are immaterial. Firstly, 
since £ zero <C D the frozen correlation length D of the 
field, and its phase, does not set the scale at which line 
zeroes appear. Equally importantly, we have a situation 
in which the density of line zeroes depends entirely on the 
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scale at which we look. Such fractal behaviour cannot be 
understood as representing vortices. We would not wish 
to identify the length scale £ zero with ^ef at this time, 
even if the latter could be defined, or D with either. 

This is not surprising. Although the field correlation 
length £(Ai) may have frozen in at D by At = 0, the 
symmetry breaking has not been effected. In a classical 
sense at least, vortices can only be identified once the 
field magnitude has grown to its equilibrium value 



(7) 



and we should not begin to count them before then. Even 
though the long-range correlation length £(At) will not 
change substantially in that time, this will not be the case 
for £ 2eT . (Ai), as long-range modes increase in amplitude 
as the field becomes ordered. 

Specifically on continuing to use (J3J) for times t > £qTq 
we see that, as the unfreezing occurs, long wavelength 
modes with k 2 < t/rQ — cq grow exponentially. Provided 
CqTq 1 they soon begin to dominate the correlation 
functions. Let 



G n (At) 



dkk 2n G(k,At) 



be the moments of G(k, At), now of the form G(k, At) = 



l TQ+At 



dt' exp 



/•e TQ+At 


fc2 +( £ o-^")) 


-2 J dt" 







We find that 



Gn(At) 



In 



2 n+l/2 



dt' 



[<' + /2/2]»+ 1 /2 : 



(8) 



where we measure the dimensionless time At in units of 
t = ^Jtq from t = €qT, q and I n = J Q dkk 2n e~ k . For 
small relative times the integrand gets a large contri- 
bution from the ultraviolet cutoff dependent lower end- 
point, and we recover As long as the endpoints make 
a significant contribution to the whole then the density 
of line zeroes derived from (j^) will be strongly depen- 
dent on scale. Only when their contribution is small and 
dn zero /dl is small in comparison to n zero /l at I = £o 
can we identify the essentially non-fractal line-zeroes with 
vortices, and £ ze ro with £de/- 

Aii and At2 , the times at which the exponential modes 
begin to dominate in the integrands of G\ and G2, can 
be determined by comparing the relative strengths of 
the contributions from the scale-independent saddlepoint 
and the endpoint in (||) . \ip\ = Aii/iandj»2 = At^/i&re 
the multiples of i at which it happens then, for the former 
to dominate the latter requires 

e^/Pi /2 > t^ /4 /V2^ and 

see that it takes longer for 
the long wavelength modes to dominate G2 than the or- 
der parameter G\ . Because of the exponential growth, pi 
and P2 are O(l). We stress that these are lower bounds. 



Since £(Ai) and ^def(At) use different attributes of the 
power spectrum G{k) e~ k 1 there is no a priori reason 
as to why they should be related. However, if the linear 
equation (||) were valid for later times At > At 2 then the 
integrals are dominated by the saddle-point at t' = At, 
to give a separation of line zeroes t; zero (At) of the form 



&™(Ai) 



Gi(At) 
G 2 (At) 



4 At 2 {tq 



1/2 



4At , , N 



approximately independent of the cutoff I for / = O(£o)- 
Only then, because of the transfer of power to long wave- 
lengths, do line zeroes become widely separated, and 
£ zero (Ai) w £de/(Ai) does begin to measure vortices. 
Thus, if the order parameter is large enough that it 
takes a long time (in units of t) for the field to pop- 
ulate its ground states (0) we would recover the result 
Cdef — 0(D) as an order of magnitude result from (g), 
but for different reasons. 

Whether we have time enough depends on the self- 
coupling [3, which determines when the linear approxima- 
tion fails. At the absolute latest, the correlation function 
must stop its exponential growth at At = At sp , when 
(|(/>| 2 }, proportional to G\, satisfies (Q). Let us suppose 
that the effect of the backreaction that stops the growth 
initially freezes in any defects. This then is our prospec- 
tive starting point for identifying and counting vortices. 

To determine At sp we need to know the strength of 
the noise. In our dimensionless units we have 

(77 Q (M)?7 6 (-k',t')) = S ab {2k B T){$/a )6(t-t')f(k-k'). 
As a result, (|7]) is satisfied when 

Gi(At v ) = n 2 altllpk B T c = ir 2 /^l-T G /T c , (10) 

and we have used the fact, for 4 He at least, that the 
absolute temperature does not change dramatically dur- 
ing the quench. Tq is the Ginzburg temperature, above 
which the potential difference between the false vacuum 
4> = and the symmetry-broken ground states is suffi- 
ciently small to be easily bridged by thermal fluctuations 
in a correlation volume. 

In order that G\ is dominated by the exponentially 
growing modes, so as to recover the modification of pre- 
diction (Q) in the form ([)]) for identifiable vortices, the 
condition At sp > Ati,At% becomes, from ( [To|) 

(t q /t q )(1-T g /T c )<Ctt\ 



(11) 



on restoring tq, where C ss 1. How strongly the inequal- 
ity should be satisfied in ([ll]) is not obvious, assuming as 
it does that the backreaction is effectively instantaneous, 
and that mean field critical indices are valid. The former 
probably means that the density above is an overestimate 
if self-consistent linearity works here as in the relativistic 
theory [HJ . On the other hand, the interaction between 
field modes will tend to redistribute the power back to 
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shorter wavelengths in the short-term, and the retarded 
nature of G(r) after At sp will impose oscillations that 
may have consequences, as well as the possible incorpo- 
ration of phonon modes. However, there is no way that 
the inequality can be remotely satisfied for He, when 
subjected to a slow mechanical quench, as in the Lan- 
caster experiment, for which tq/tq = O(10 10 ), since the 
Ginzburg regime is so large that (1 — Tg/T c ) = O(l). 
That is, field growth must stop long before vortices are 
well defined. Further, since the experiment leaves the 
superfluid in the Ginzburg regime, thermal fluctuations 
will inhibit the creation of stable defects, although it may 
be that incoherent fluctuations, even if not vortices, will 
give a signal. [In no experimemt is the vortex separation 
given by the corresponding Ginzburg correlation length 
£g = 6o/\/l — Tq/T c , independent of tq, as was orig- 
inally thought jUl]. Adopting renormalisation improved 
critical indices cannot repair such a deficit. Whatever, 
there is no reason to expect a vortex density (||) . 

The situation for 3 He — B is potentially very different, 
with a much more complicated 50(3) x SO(3) x £7(1) 
order parameter A a i |[6[ permitting several types of vor- 
tex, whose cores are not necessarily normal. However, 
for rapid quenches a TDLG approach is valid fl7| with 
potential V(A) = a(t)\A al \ 2 + 0(A A ). Prior to back- 
reaction setting in, this is similar to the He case consid- 
ered above in tmj , although self-consistent linearity would 
be somewhat different. Suppose that (|TT|), which looks 
like an inequality between the quench rate and the equi- 
libriation rate that permits unstable modes time to dom- 
inate over short-range fluctuations, remains true, albeit 
with new coefficients. Firstly, for 3 He — B, the Ginzburg 
regime is very small, with 1 — Tq/T c = O(10~ 8 ) and 
mean-field is a good approximation. Secondly, in gen- 
erating the phase transition by nuclear reactions 
rather than by mechanical expansion, the quench rate is 
increased dramatically, with tq = 10 2 — 10 4 . The in- 
equality ( pi] ) is satisfied by a huge margin and we can 
understand the success of the Helsinki and Grenoble ex- 
periments. We note that, if the density is given by (j|), 
then this overestimates the primordial density (H) as 



ndef 



(12) 



where p = At sp /t > 1 = 0(1). 

Much more detailed numerical modelling in three spa- 
tial dimensions is needed before we can draw any quan- 
titative conclusions at later times. In d < 3 dimensions 
Gi and G2 are replaced by Gu-X)/2 and Gu+i)/i- It 
is easier to create defects in less than three dimensions 
(particularly in one dimension jig] ) since, with smaller 
n, the saddle point in (JsJ) dominates more easily over the 
weakened ultraviolet endpoint. 

However, we believe that we have presented a convinc- 
ing case that, when line-zero density is well-defined, it 
can be equated to vortex density. Because line-zero den- 
sity is based on different attributes of the power G(k) 



from the effective long-range correlation length £(Ai), 
we should not equate them directly. Only if there is time 
enough before the field populates the ground states do 
we recover the result (|l2|) , a specific case of (||) , although 
for different reasons. Moreover, the necessary condition 
([ll]) discriminates cleanly between the experiments on 
4 He and 3 He. More details will be given elsewhere |J. 

Finally, a constraint like ([n]) is not specific to the 
TDLG theory. A similar approach to weakly coupled 
relativistic quantum field theory |l5[ permits a compari- 
son with Kibble's predictions After some rearrange- 
ment, it can be seen Jj"5| that a necessary condition 
for the appearance of well-defined 0(2) vortices after 
the immediate implementation of the transition is that 
(tq/to) 2 (1 — Tg/T c ) < 1, although we did not appreciate 
this at the time. 
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